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Bi-spectral satellite thresholds for precipitation spec- 
ification are explored with visual and infrared satellite 
data collocated with Service-A hourly observations for 137 
surface stations in the southeastern tJnited States. The 
data span the month of August 1979 and total 70,623 observa- 
tions, including 538 daylight precipitation observations. 

The distributional and statistical differences of four 
satellite resolution sizes ranging from 484 to 2025 nmis are 
explored and determined to be significant in the representa- 
tion of weather conditions. Precipitation and no-precipita- 
tion data can be statistically differentiated with the 
visual and infrared mean and standard deviation values. 

For overcast ceiling reports, a simple linear bi-spec- 
tral threshold based on a 50% probability of precipitation 
is defined as extending from albedo 1.00 to 0.60 with asso- 
ciated cloud top temperatures 290K and 210K, respectively. 

For overcast and broken ceiling reports, an albedo greater 
than 0.80 specifies a 50% probability of precipitation. 
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INTRODUCTION 



The determination of precipitation occurrence and 
amounts is an important factor in scientific, commercial, 
and operational endeavors. Scientific uses are concentrated 
in the fields cf meteorology, hydrology, and oceanography, 
where precipitation is essential in analysis, diagnosis, 
prediction, and verification. within meteorology, precipi- 
tation serves as both a forcing and response element in the 
study of daily weather and climatology. Indeed, precipita- 
tion is a critical input for climate research and into gen- 
eral circulation models which premise to extend the time 
frame of skillful weather forecasts. Commercial uses encom- 
pass agriculture, forestry, transportation, communications, 
water resource management, and many others. 

Despite the importance of precipitation data to a vari- 
ety of fields, there are serious shortcomings in current 
precipitation determination. These shortcomings are due to 
areal and economic limitations imposed upon the land-based 
rainfall monitoring systems. k possible solution is embod- 
ied in precipitation information extracted from satellite 
data. With the advent cf high resolution, multi-spectral 
channel satellites in the late 1970's, satellite derived 
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precipitation data are being studied as a viable method zo 
complement and supplement conventional rainfall data. 

The satellite image interpreter does a subjective analy- 
sis based on the gray shade variations, representing a range 
of digital counts, that appear in the satellite image. How- 
ever, satellite data contain more information within the 
digital values than can be resolved by the human eye in pho- 
tographic images. The satellite digital counts input into a 
computer allow use of the full range of the digital values. 

Until recently, computer processing of satellite data 
has been confined largely to research uses. Acquisition, 
storage, and processing of the huge volumes of digital sat- 
ellite data could not be handled operationally in real time. 
However, with the recent advent of more capable mini-com- 
puter systems, such as the United States Navy's Satellite 
Data Processing And Display System (SPADS) developed by the 
Naval Environmental Prediction Research Facility (NEPRF) 
Monterey, California, real time quantitative use of digital 
satellite data has become a reality. With the operational 
availability of such systems as the SPADS unit, there is a 
need for numerical schemes to aid in the objective 
specification of current weather conditions. 
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This thesis concentrates on the specification of visual 
and infrared satellite data thresholds in determining pre- 
cipitarion occurrence and qualitative precipitation rates in 
a mid-latitude coastal environment. The data set used con- 
sists of collocated Geostationary Operational Environmental 
Satellite--East (GOES-E) satellite data and hourly surface 
observations at East Coast and Gulf Coast United States sta- 
tions, south of 40®N, for the month of August 1979. 

The use of satellite data for precipitation specifica- 
tion is not new. There is the recognized limitation that 
infrared and visual satellite sensors are measuring proper- 
ties associated with small cloud particles and not precipi- 
tation sized particles. Nonetheless, Muench and Keegan 
(1979) specified quantitative precipitation rates, Liljas 
(1981a, 1981b) specified qualitative precipitation rates, 

and Love joy and Austin (1979) delineated rain versus no-rain 
cases using visual and infrared satellite data. Del Beato 
(1981) used cloud top temperatures within a restricted cloud 
case classification to derive qualitative precipitation 
rates. 

The currently available precipitation study results are 
based on data sets with region, season, and size 
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li mitarions. 



This research effort will use data from 



stations covering more than 420,000 square nautical miles 
(nmi2) in the eastern and central United States with a total 
of 70,623 observations (538 precipitation observations) . In 
comparison, the relatively comprehensive precipitation study 
of Muench and Keegan (1979) was based on 552 cases (300 
rainfall cases) from five stations in the northeastern 
United States for April through November 1977. The signifi- 
cantly larger size of the present sample will allow better 
statistical determination of appropriate distributions and 
threshold value significance. 

The primary objective of this thesis is to investigate 
specification of precipitation versus no- precipitation from 
satellite visible and infrared digital counts. Addition- 
ally, in precipitation cases, the feasibility of qualitative 
specification of light versus moderate/heavy precipitation 
and guantiative specification of convective versus continu- 
ous precipitation are investigated. 

The thesis is organized into five chapters. Chapter II 
reviews the satellite data based precipitation studies. 
Chapter III describes the data set, the data processing and 
the testing program. Chapter IV describes the results. 
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chapter 



V states the 



conclusions and 



suggest s 



further 



research . 
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II . PRECIP ITATION SPECIFIC ATION 



A. INTRODUCTION 

The specification of precipitation and tha estimation of 
rainfall rates using satellite imagery have bean studied 
using a variety of methods over a wide spectrum of time 
scales. This review will concentrate on those methods 
developed for synoptic scale and mesoscale analysis of pre- 
cipitation on a diurnal or shorter time scale. The methods 
reviewed include bi-spectral and infrared threshold (Muench 
and Keegan, 1979; Liljas, 1981a, 1981b; Lovejoy and Aus- 
tin, 1979; Del Beato, 1981; Hylie, 1982) and life history 
(Scofield, 1981; Griffith et al. , 1978; Stout et al. , 

1979; Wylie, 1979; Negri and Adler, 1981). 

B. BI-SPECTRAL AND INFRARED THRESHOLD 

The bi-spectral threshold method, in which infrared and 
visual satellite data are used, involves mapping the extent 
and distribution of precipitation. Combining the visual and 
infrared data provides information on the cloud temperatures 
(infrared data) and on the cloud thickness (visual data). 
Thus, while use of the visual or infrared data alone may 
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the combina- 



have limitations in specifying precipitation, 
tion of both may succeed at specifying precipitation. The 
multi-spectral satellite channels introduced on satellites 
in the late 1970*s yielded the possibility of bi-spectral 
thresholds. Threshold values and study condition parameters 
of selected bi-spectral studies are summarized in Table I. 

TABLE I 

Summary of Bi-spectral and Infrared Threshold Values 



STUDY 


LOCATION 


asEs 


TIME OF 
YEAR 


THRESHOLD 

INFRARED 


VALUE 

VISUAL 


Muench and 
Yeegan 
(1979) 


Northeastern 
United States 


552 

obser- 

vations 


April- 

November 

1977 


-12®C 


0,60 


Liljas 

(1981) 


Scandinavia 


- 


May 1979 
August 
1979 


-12®C to 
-15®C 


- 


Lovejoy and 


Montreal 


17 days 


June 1977 


-21®C, 


.80, .88 



Austin -26^C, 

(1979) -41®C 

* Visual threshold based on normalized scale from 0-1 

Muench and Keegan (1979) studied precipitation specifi- 
cation using GOES visual and infrared satellite data and 
hourly rainfall climatological data for five stations in the 
northeastern United States for the period April through 
November 1977, Their data set consisted of 552 
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observations, comprised of 300 rainfall observations and 252 
cases of either nonprecipitating cloudy or fair weather 
observations. The visual (1 km resolution) and infrared (7 
km resolution) satellite data were area averaged over 7x7 
square kilometers (km^) and 14 x 14 km^, respectively. A 65 
point visual data array (8x8 plus the cenrer point) and a 
17 point infrared data array (4x4 plus the center point) 
were centered ever each station. The GOES visual data were 
normalized using reflection values from Liou (1976) with the 
modification of lower absorption and higher transmission to 
compensate for Liou’s treatment of the complete solar spec- 
trum. The anisotropic radiation of clouds was corrected 
with functions calculated by Mue nch and Keegan (1979) from 
ground-based radiometers and satellite measurements. From 
these data, they determined probabilities for precipitation 
greater than .01 and .10 inches in one hour and the amount 
of precipitation for the hour following the satellite 
observation (see Figs. 1, 2, and 3). 

Muench and Keegan (1979) did not provide the standard 
deviations for the data in these figures. However, they 
stated there was ’’considerable uncertainty in the specifica- 
tion of rainfall amount." As an example, they stated that 
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IR Cloud Tcmptrofurt (DopC) 

Figure 1. One-hour Rainfall as a Function of Normalized 

Cloud Reflectivity and Infrared Cloud Temperature 
(from Muench and Keegan, 1979) 




Figure 2. Probability of One-hour Rainfall Greater than or 
Equal to 0.01 inches as a Function of Cloud 
Reflectivity and Infrared Cloud Temperature (from 
Muench and Keegan, 1979) 
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1.0 




Figure 3. Probability of One-hour Rainfall Greater than 

0. 10 Inches as a Function of Cloud Reflectivity 
and Infrared Cloud Temperature (from Muench ana 
Keegan, 1979) 

using Fig. 3 "for a one-hcur rainfall specification of 0.10, 
two- thirds of the values would fall between 0.25 and 0.04." 
Muench and Keegan stated that their figures emphasize the 
requirement fcr both visual and infrared data to specify 
precipitation amounts. 

Liljas (1981a, 1981b) developed a bi-spectral cloud 

classification based on visual and infrared data from the 
polar orbiting TIR0S--6 satellite (see Fig. 4 and Table II). 
The data set consisted of a limited number of daily observa- 
tions, chosen for their synoptic characteristics, in May and 
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Ch4 




Figure 4. Two Dimension Decision Space for Typing Clouds 

from Visual and Infrared Digital Counts (Table II 
defines the symbols for the clouds) (from Lilias, 
1981a) 



August 1979 over a region encompassing Norway, Sweden, Fin- 
land, and the Baltic Sea with weather charts providing the 
ground truth. Based upon the precipitation threshold 
results of Muench and Keegan (1979), Liljas chose a cloud 
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TABLE II 



Cloud Classification to be used with Fig, 4 
(from Liljasr 1981a) 



Main and Cloud typaa: 
1« Cuzauloniabus 

2, Hlabostracua 

3, Cirrostracus 



4 « Cumtilus congestus 



5» Stracocumulus 



6. Haze/Scracua 

7. Land 



8« Water 



a 




f 




1 

g 

1 

J 

1 , 



k 



Stora cloud with high top 

Squall cloud with scattered showers 

Large vertical chicicnesa 
Rather low topside 

Dense clrrostratus 
Cirrus 

Thin cirrus over water 
Dense altostracus 
Thin altostracus over water 
Thin altostracus 

Dense altocumulus 
Large piled up cumulus 
Rather snail and flat conulus 

Dense stracocumulus 
Ordinary 

Slightly piled up c^xmulus 
with clear areas in between 
Very dense haze/scracus 
Dense haze/scracus 

Ordinary haze/scratus 
Cumulis humilis 



Haze over water 

Planting season spring or autumn 
Warm green season 

Cold 

Vara 
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top temperature threshold of -,12®C to -15®C to classify 
cumulonimtus and nimbostratus clouds. Starting with this 
cloud classification and the assumption that the highest and 
densest clouds produce the maximum precipitation amount, 
Liljas suggested a qualitative precipitation intensity scale 
based on the sum of the visual and infrared satellite digi- 
tal counts (see Table III) . These sums represent the areas 
of the Liljas nimbostratus and cumulonimbus cloud types in 
his bi-spectral cloud classification (see Fig. 5). 

TABLE III 

Threshold Values Describing Precipitation Intensity Levels 
as Applied in Fig. 5 (from Liljas, 1981a) 

T he S um of Di gital Levels 

Ch 1 + Ch 4 391-310 light rain 

31 1-330 
331-350 
35 1-370 
37 1-390 

390 very strong rain 

Lcvejoy and Austin (1979) studied rain mapping of cloud 
areas based on GOES visual and infrared satellite data over 
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Figure 5. 



Precipitation Intensity Classification from 
Visual and Infrared Digital counts. The 
Precipitation Area is Represented with the Dark 
Diagonal Lines. (See Table III for the 
mathematical description of the intensity areas.) 
(from Liljas, 1981af 



Montreal, Canada, and the trop 
pheric Research Program Atlantic 
data) with radar data providing 
real data set consisted of 17 o 



ical Atlantic (Global Atmos- 
Tropical Experiment, GATE, 
the ground truth. The Mont- 
bservations over three days 
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during June 1977. Working with 4 x 4 km resolution satel- 
lite image. Love joy and Austin plotted two dimensional fre- 
quency grids for the radar-determined rain and no-rain 
points on a 25 X 25 array (see Figs. 6 and 7). The visual 



data were normalized by selectin 
mest” values in each image and 
radiances between 0 and 1. 

Lovejoy and Austin (1979) s 
cumulus rain data distribution o 
bution was to a good approximat 
sian,” They do not describe o 
support this assertion. The no- 
were described as a bimodal dis 
the low visual and low infrared 
near the rain peak but shifted s 
In most cases, the separation 
rain cases was statistically si 
ity ranging from to 50% 

samples came from the same popul 

The Lovejoy and Austin (1979 
plots for non-cumulus storms wer 



g the "brightest” and "dim- 
linearly interpolating the 

tate, with reference to the 
f Fig. 6 that, "The distri- 
ion a two-dimensional Gaus- 
r provide the statistics to 
rain cumulus cases (Fig. 7) 
tribution with one peak near 
values and the other peak 
lightly toward lower values, 
of the oumulus rain and no- 
gnificant with the probabil- 
that the rain and no-rain 
ation. 

) two dimensional frequency 
e limited ~o one case. The 



significant differences between the cumulus 



and non-cumulus 
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Figure 6. Frequency Plot of Rain Distribution for GATE day 
248^ 1300 GMT (from Lovejoy and Austin^ 1979) 
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Figure 7. Frequency Plot of No-rain Distribution for GATE 
day 248, 1300 GMT (from Lovejoy and Austin, 1979) 
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data sets were that the non -cumulus no-rain plot lost its 
bimodal character, relative to the cumulus no-rain plot, and 
appeared as a broad two dimensional Gaussian distribution. 
The non-cumulus rain plot points fell within the no-rain 
distribution, but were shifted slightly higher in the vis- 
ual. The separation of the non-cumulus rain and no-rain 
cases was not statistically significant, with greater than a 
50% probability of the rain and no-rain samples coming from 
the same population. 

Lovejoy and Austin (1S79) attempted to further classify 
the cumulus rain and no-rain cases into no-rain, light rain, 
and heavy rain. Rainfall rates greater than 2 mm-h-^, as 
determined by radar, were defined as heavy rain. As 
expected, the mean of the heavy rain cases was shifted 

slightly towards higher visual and infrared values than the 
mean of the light rain cases. However, rhe shift was so 
small that there was at least an 80% probability of the 
light rain and heavy rain cases coming from the same popula- 
tion. Lovejoy and Austin (1979) concluded that "little if 
any rainfall-rate information is contained in a single (vis- 
ual and infrared) satellite image." 
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Lovejcy ar.d Austin (1979) tested a spectral threshold 
technique for rain area mapping. Each satellite image of 
400 X 400 km was divided into one hundred 40 x 40 km boxes. 
The 100 subareas were each checked with radar to determine 
the total number of rain areas. An equal rotal number of 
satellite subareas were classified as rain areas. The sat- 
ellite subareas with the highest visual and highest (cold) 
infrared values ware classified as rain areas, until the 
total number cf satellite rain areas equaled the total num- 
ber of radar determined rain areas. This spectral threshold 
technique was applied to three days accumulation of data and 
is shown in Tables IV and V. When compared with the success 
of the two dimensional frequency plot method, the visible 
and infrared thresholds averaged 45% and 58% worse, respec- 
tively. The accuracy of the visual threshold is limited by 
the extent of low, thick clouds and the infrared threshold 
is limited by the extent of the cirrus clouds in the satel- 
lite image. Lovejoy and Austin (1979) concluded that ”the 
errors involved in using a 'best threshold* are very large 
indeed. " 

Del Beato (1981) studied correlations between cloud top 
temperatures (based on NOAA-5 satellite data) and rainfall 
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TABLE IV 

Statistical Comparison of Rain Area Mapping Techniques 

(R /R X 100 indicates "percentage of correct satellite rain") 
(from Lovejoy and Austin, 1979) 





OpC. 2*D Boundafy 


IR Optimum Threshold 


Visible Optiimim Threshold 


Area Rain 
Coverafe CO 


ToUlNo. 
of Point! ^ 


Day 


X I00» 


IR(K) 


(Ma/M) X 100 


(Scale: 0-1) 


(iVa/R) X 100 


c^r€ 


242, 24J. 246 
247, 248. Ul 
252. 261 


65 


<232 


33 


>0.68 


64 


13.8 


47706 




152 


36 


<232 


20 


>0.88 


31 


9.7 


. 40361 


110 


36 


<247 


33 


>0.80 


48 


24.0 


33738 


133 ‘ 


33 


<234 


32 


>0.81 


49 


13.9 


22331 



TABLE V 

Statistical Comparison of the Accuracy of Rain Areas 
(from Lcvejcy and Austin, 1979) 









Number of Images 




Error 






Technique 


Region 


or Sequences 


Bias 


Factor 






2-D Pallern 


Montreal 


17 


1.13 


1.26 


0.22 




Matching 
2-D Pattern 


Montreal 


3 


1.08 


1.19 


0.18 




Matching 
Optimum IR 


Montreal 


3 


1.38 


1.74 


0.71 




Threshold 
Optimum Visible 


Montreal 


3 


1.54 


1.59 


0.58 




Threshold 
2-D Pattern 


GATE 


8 


1.21 


1.41 


0.25 




Matching 












totals for 30- and 


60 -min intervals over 


eastern 


Australia . 


The 


satellite data 


had a 


60 km^ maximum 


rs 


solution 


at subsa- 


tellite point and 


cloud 


top temperatur 


93 


were area averaged 


fo r 


a resolution 


of 200 


km2. The 21 


dara sets 


were first 



classified according to synoptic situation in a rough 
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attempt, to group the data by cloud type, droplet spectra, 
and air mass trajectory. The initial results suggested that 
the cloud top temperature determined an upper limit on rain- 
fall amount, with the maximum increasing as the cloud top 
temperature decreases. A linear correlation analysis to 
determine a quantitative relationship between rainfall 
amount and cloud top temperature gave indefinite results. 

Further study of surface and radiosonde observations 
indicated that classification by proportion of cumuliform 
cloud reports to all cloud reports and subcloud layer humid- 
ity might be mere appropriate (Del Beato, 1981). This clas- 
sification resulted in a correlation coefficient of 0.90, 
excluding cases with cumuliform portions less than 50% and 
dew-point depressions of greater than 6®C. Finally, a com- 
posite frequency distribution was calculated based on three 
cases, all southwesterly stream situations described as 
"post-frontal cellular convection cases in cyclonically 
curved flow." The fitted equation was: 

f = 0. C57 - 0.004CTT - 0.054R (1) 

where f is the rainfall frequency, R is the 30-min rain 
total (mm) , and CTT is the cloud top temperature (°C) . The 



32 



equation was fitted to 41 independent f values. This equa- 
tion is associated with a correlation coefficient of 0.79 at 
the 99% confidence level. Equation (1) indicates no rain 
from clouds warmer than +13®C and a maximum 30-min rainfall 
of 2.5 mm for a cloud top temperature of -20®C. 

In summary, Del Beato (1981) found than cloud top temp- 
eratures and 30- and 60-min rainfall torals indicated sta- 
tistically significant relationships for cloud systems with 
a high proportion of cumulus clouds and high subcloud humid- 
ity. Additionally, as cloud top temperatures decrease to at 
least -35®C, rainfall totals increase. 

Wylie (1982) attempted to correlate rainfall occurrence 
with radiosonde soundings, hourly Service-A observations, 
and visual and infrared satellite data. His data sample was 
restricted to '’large-scale cloud cover" areas with wide- 
spread precipitation (rain gauge reports varied less than 
20%) for the Great Plains States region for the period 27 
February 1981 through 4 January 1982. From thirteen parame- 
ters derived from the three data sources (see Table VI) , the 
best linear regression equation for estimating rainfall 
ra te s wa s : 

6-hour rain (in) = 1.0242 + 0.380Pw - 0.0304Qc 

- 0. 0047Ct (2) 
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where Pw is the vertically in 
vapor (in ) , Qc is the moisture 
Ct is the cloud rop temperature 
a linear correlation coefficient 
equations were also determined 
alone and for a combination of P 
all three data types were availa 
regression equation was: 

6 hour rain (in) = 2.10-0 

The correlation coefficient was 
that the synoptic scale data b 
suited for estimating broad chan 
ciated with changes in air mass 
mating rainfall rates associate 
processes . 

C. LIFE HISTOEY 

The life history methods are 
itation estimation schemes ba 
first, that significant rainf 
clouds, and second, that convect 
and measured in satellite imag 
manual analyses of convective cl 



tegrated precipitable water 
convergence (g/kg/day) , and 
(Kelvins) . Equation (2) has 

of 0.60. Linear regression 
for the three parameters 
w and Qc to be used when not 
ble. The cloud temperature 

.OOSCt (3) 

-0.35. Wylie (1982) stated 
ase measurements were best 
ges in rainfall rates asso- 
es and not suited for esti- 
d with small scale dynamic 

empirically derived precip- 
sed upon two assumptions, 
all comes from convective 
ive clouds can be identified 
es. These methods involve 
oud areas in a sequence of 
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TABLE VI 

Correlaticn Coefficients for Determination of Precipitation 
Based on the Three Data Types (from Wylie, 1982) 



NEASORED 

PARAMETI^ 


CORR ELATION 
WITH 6 HOUR 
PRECIP. REPORT 


NUMBER 
OF CA^S 


F 


Vertically integrated 
precipitacle water vapor 




0.48 


196 


58* 


Cloud top brightness 




-0.44 


184 


44* 


Cloud top height 




-0.40 


190 


36* 


Moisture convergence 




0.38 


184 


3 1* 


Cloud top temperature 




-0.35 


199 


27* 


Bubble model predicted cond. 




0.27 


115 


9 


500 mb vorticity advection 




-0.21 


173 


8* 


Parcel lifted index 




-0.20 


200 


8* 


700 mb temperature advection 




0.20 


173 


7* 


Sfc temperature advection 




0.19 


156 


6 


850 mb temperature advection 




0.17 


189 


6 


Wind convergence (sfc) 




0.09 


167 


1 


Vertical wind shear 




0.03 


156 


0 


* Significant correlation at the 


! 9955 lave 


1. 




visual, infrared, or both 


visual and 


infrared sat 


ellite 


images. Threshold values 


and 


study condition parameters 
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associated with published life history studies are summa 
rized in Table VII. 



TABLE VII 

Summary of Life History Threshold Values 



STUDY 


LOCATION 


CASES 


TIME OF 
YEAR 


THRESHOLD 

INFRARED 


VALUE 

VISUAL 


Griffith, 
et al. , 
(1978) 


Florida, 

Venezuela, 

Honduras, 

and hurricanes 

impacting East 

Coast United 

States 


34 days 


summers 

1969-1976 


-20°C 


80 counts* 


Stout, 
et al., 
(1979) 


tropical North 
Atlantic 


57 obser- 
vations 


September 

1974 


-26°C 


0.45 albedo 
(sun over- 
head) 


Wylie 

(1979) 


Montreal 


6 days 


June 1977 
September 
1977 


-16®C 


- 


Negri and 
Adler 
(1981) 


Oklahoma , 
Arkansas, 
Missouri 


1 day 

(15 thunder- 
storms) 


April 24, 
1975 


-27®C 


- 



* ATS-3 sateUite 



The Scofield/Oliver (Scofield, 1981> analysis follows a 
decision tree procedure to estimate half-hourly rainfall for 
deep convective systems within tropical air masses. Using 
enhanced infrared and high resolution visual satellite 
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images, the technique involves first identifying the active 
convective portion of the cloud, or cluster, from two con- 
secutive satellite images. Once the acrive portion is iden- 
tified, the half-hourly rainfall estimation is computed 
based on such factors as cloud top temperature, cloud 
growth, and departure of precipitable water from a summer- 
time normal. 

The Grif fith/Hoodley (Griffith et al. , 1978) technique 

is designed to estimate rainfall in the tropics, over large 
space and time scales, using geosynchronous visual or infra- 
red satellite imagery. This time-dependent technique was 
empirically derived as a relationship between cloud area, 
echo area, and rain rate for two areas in south Florida, 
with raingage-radar providing the ground truth, and was then 
tested in other tropical areas. This scheme was subse- 
quently tested further in extratropical areas (Griffith et 
al., 1980), with modifications to the rainfall amount 
predicted. 

The determination of a cloud area-rainfall relationship 
first required the specification of both a visual and an 
infrared threshold to define the cloud area. The visual 
brightness threshold, normalized for radiation geometry, was 
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80 counts for the third Application Technology Satellite 
(ATS-3) and the infrared threshold was 253K (-20°C) . The 

thresholds were based on a comparison of the clouds with a 
given maximum digital count and the radar echoes associated 
with these clouds. 

The empirical cloud a rea -rainfall relationship was 
derived as a two step process. First, a relationship 
between the cloud area and the radar echo area, normalized 
for the maximum area achieved by the cloud or cluster, was 
established for the visible and infrared satellite data. 
Second, the relationship between the echo area and rain vol- 
ume was determined and was of the form: 

Ev = I Ae (5) 

where Rv is rain volume per hour (m^-h” ^) , I is rain in units 
of ( m3-km~2-h“ 1) , and Ae is the echo area (km2) defined by 
the 1 mm-h“i rain rate. Thus, given a rime sequence of con- 
vective clcuds (or cluster areas) measured from visible or 
infrared satellite images, volumetric rain rate can be 
estimated . 

Stout et al. (1979) modified the Griffith/Woodley tech- 
nique (Griffith et al., 1978) ro estimate volumetric rain 
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rate directly from a cumulonimbus cloud area and area change 
according to the equation: 

B = a^A + ai dA/dt (6) 

where R is the volumetric rainfall of the cloud (m^-s-i), A 
is the cloud area (mz) , dA/dt is the change of cloud area 
over time (mz-s~‘), and aQ and a, are constants with dimen- 
sions m-s-i and m, respectively. The two constants were 
calculated by a least squares fit of cloud area-rain rate 
pairs based on visible and infrared geosynchronous satellite 
data and 5.3 cm ship radar rain data collected during GATE. 
The cloud area and its change are defined by the threshold 
value. The visible threshold for cloud area calculations 
was 60 digital counts on the ATS-3 (corresponding to an 
albedo of 0.45 with the sun overhead), or equivalently 172 
digital counts on the first Geosynchronous Meteorological 
Satellite (SMS 1). The infrared threshold was 160 digital 
counts (-26°C). The standard error between the estimated 
rainfall and the mean radar rainfall was 62% and 76% for the 
visual and infrared equations respectively. 

Wylie (1979) attempted to use the tropical convective 
rainfall techniques of Griffith et al. , (1978) and stout £t 
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al., (1979) for estimating precipitation in Montreal, Can- 
ada. Using visual satellite data, corrected for the chang- 
ing sun angle (Mosher, 1975), infrared satellire data, and 
10.0 cm radar meaured rainfall rates, Wylie studied six days 
of precipitation, three days each in June and September 
1977. Wylie concluded that because of air mass differences 
between Montreal and the tropics, the Griffith and Stout 
estimation techniques did poorly in Montreal, Canada. The 
singlemost important limitation with these two schemes was 
the difficulty of measuring cumulonimbus cloud area when the 
"anvils were often merged into large cloud masses and the 
extensive stratus cloud cover often obscured the pictures." 
Wylie also noted that the Griffith el. al. (1978) threshold 
of -26®C had to be changed to -1 6°C for the summertime Mont- 
real, Canada, area. With the warmer cloud top temperatures 
the cloud areas were a larger, more appropriate size for 
tr acking. 

Wylie (1979) then attempted to combine sounding data 
input into a one- di mens iona 1 model (Simpson and wiggert, 
1969) and satellite cloud cover measurements to estimate 
rainfall for Montreal. With the GATE measurements for rain 
rates associated with sat ellite- derived cloud areas and the 
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model output, rainfall rates were eszimated by multiplying 
the two values. The most accurate estimations were for the 
cumulus clouds in the warm air masses occurring in June, the 
cases the model was designed to handle. Wylie concluded 
that in order to estimate rainfall in all geographical areas 
and seasons a more sophisticated model would be needed. 

Negri and Adler (198 1) did one case study of fifteen 
thunderstorms in the Oklahcma, Arkansas, and Missouri area 
on 24 April 1975. They used radar data for ground truth and 
had special 5 minute GOES-E satellite passes over the area 
of interest. They were able to determine that the precipi- 
tation began falling, as indicated by radar data, for cloud 
zop temperatures ranging from 229K to 260K {-44®C to -13°C). 
The mean cloud tcp temperature value was 247K (-26°C) . 
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Ill 



CAT A PR OCESS ING 



A. INTRODUCTICN 

The data set assembled for this study consists of collo- 
cated GOES-E satellite data and Service-A hourly surface 
observations for the southeastern United States during 
August 1979. The GOES-E data consists of 10 x 10 pixel 
matrices of visual and infrared satellite dara centered over 
each of 137 surface stations (Fig. 8) all south of 40^N. 
The satellite data are measured with the Visual Infrared 
Spin Scanned Radiometer (VISSR) which have subsatellite 
point spatial resolutions of 1 and 7 km for the visual and 
infrared channels, respectively. The GOES-E navigation was 
completed by Han-computer Interactive Data Access System 
(McIEAS) at the University of Wisconsin using the full reso- 
lution visual data, with an accuracy of 1-2 pixels (1-2 km). 
The full resolution visual data were averaged to a 7 km res- 
olution, to equal the infrared data resolution. The visual 
and infrared digital counts range from values of 0-255. The 
10 X 10 pixel GOES-E visual and infrared satellite data each 
cover an area 45 nmi x 45 nmi at 30<*N (60 nmi x 60 nmi at 
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420N) . 



The 



Service-A hourly reports total 70,623 
observations. No Service-A specials or record-specials are 
incl uded. 




Figure 8. Geographical Locations of Service-A Station 
Report Data 



The data are divided into two no-precipi~ation catego- 
ries (Table 7III) and seven precipitation categories (Table 
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IX) to investigate pracifitation specification, convective 
versus continuous precipitation specification, and qualita- 
tive specification of light versus moderate/heavy precipita- 
tion. The pixel array size is also varied from the 10 x 10 
array size to an 8 x 8, a 6 x 6, and a 4 x 4 array size to 
investigate the differences in the data resulting from vari- 
ous resolution sizes within a particular weather condition 
cl assification . 

B. DATA SORT 

For the combined visual and infrared threshold specifi- 
cation of precipitation, satellite data for 1200-2000 GMT, 
corresponding to 0800-1600 EDT, were sorted into precipita- 
tion and no-precipitation groups (Fig. 9). The 0800-1600 
EDT interval was chosen tc avoid distortion of the visual 
satellite data due to a low solar elevation angle. The vis- 
ual data were normalized and converted to albedos based on 
the work done by Muench and Keegan (1979). This scheme cor- 
rects for the varying zenith angle as well as adjusting the 
visual satellite data for anisotropic scattering as related 
to the zenith angle. (See Appendix A for further specific 
information ccnceming the Muench and Keegan normalization 
sche me. ) 



44 




Figure 9. Flew Chart of Data Processing 
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While albedo values cannot exceed 1.00, the Muench and 



Keegan (1979) scheme allows the values to overshoot 1.00, up 
to a value of 1.20. Therefore, the visual satellite values 
are not true albedos, but estimated albedos. The extended 
visual normalized data scale was used to facilitate compari- 
son of the results in this effort to the most extensive bi- 
spectral threshold precipitation specification of Muench and 
Keegan (1979). The Muench and Keegan (1979) normalization 
scheme specifies that any computed albedo greater than 1.20 
be set equal to 1.20 to limit the unreasonably large values. 
Similiarly, the scheme specifies computed albedos less than 
0.15 be interpreted as the ground or water surface reflec- 
tance and the value 0.00 be assigned. The infrared data 
were processed in digital coun-s and converted to cloud top 
temperatures prior to statistical computations and graphical 
displays. 

The no-pr scipitat ion data (Table VIII) are comprised of 

the digital visual and infrared 10 x 10 pixel arrays of 
those Service-A station reports not showing any ”R" in the 
current weather group. Thus, the no-precipitation group 
includes stations reporting drizzle (weather codes L-, L, 
and L+) . 
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TABLE 7 III 

Classification of Nc-Prscipitation Data Groups 



Service- A 



Categor y 



Name 



Current Weather/ 
Cloud Group 



Re por ts 



2A 



No-Precipitat io n. 
Overcast Ceiling 



no ”R''/ 
Cloud Groups 
300, 030, 003 
130, 230, 013 



r 

f 

9 



1976 



023, 103, 203 
113, 123, 213 
223 



2B 



No- Precipitation, 
Overcast and Broken 



no ''R''/ 
Cloud Groups 
(above groups) 




7358 



102 , 202 , 112 , 

122 , 212 , 222 



The no-precipitation data are divided into two catego- 
ries, overcast ceiling (category 2A) and overcast and broken 



digit cloud group in the Service-A surface observation. The 
first digit indicates the amount of low clouds, where 0 is 
defined as clear, 1 is scattered (one-eighth to four- 
eighths cloud cover), 2 is broken (five-eighths to seven- 

eighths cloud cover) , and 3 is overcast (eight-eighths cloud 
cover). The second and third digit indicate the amount of 
middle and high clouds, respectively. The same 0-3 values 
defined for the low clouds are used for middle and high 
cloud amount. 



ceiling (category 2B) 



Cloud cover is based on the three 
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The 


precipitation data (Table IX) are comprised 


of the 


visual 


and infrared 10 x 


10 pixel arrays of those Service-A 


station 


reports showing 


any " R" in the current 


weather 


group . 


Two precipitation observations were excluded from 


the dat 


a set because each 


report also indicated clear 


skies . 






TABLE IX 






Classification of 


Precipitation Data Groups 




Category Name 


Servioe- A 
Current Weather 


Reports 


1 


Precipitation 


any "R" 


538 


1 A 


Precipitation any "R” and overcast 

Overcast Ceiling ceiling (as defined in 

Table VIII category 2A) 


329 


IB 


Precipitation 
Overcast and 
Broken Ceiling 


any "R" and overcast 
and broken ceiling (as 
defined in Table vIII 
category 2B) 


534 


1C 


Continuou £ 
Precipitation 


R~ f R , R+ 


1 12 


ID 


Convective 

Precipitation 


RW-, RW, RW+,TRW-, 
TRW, TRW+, TR-, TR, TR 


426 


1 E 


Light 

Precipitation 


R-, RW-, TRW-, TR- 


464 


1 F 


Moderate/H eavy 
Precipitat ion 


R, R+, RW, RW+, TRW, 
TRW+, TR, TR 


74 



The general precipitation data (category 1) are divided 
into six groups: precipitation overcast ceiling (category 
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IA) , precipitat ioa overcast and broken ceiling (category 

IB) , continuous (category 1C), convective (category ID), 
light (category IE), and modera te/hea vy (category IF) pre- 
cipitation. These seven precipitation groups are used to 
investigate precipitation specification, convective versus 
continuous precipitation specification, and qualitative 
specification of light versus mo derat e/heavy precipitation. 

C. STATISTICAL TREATMENT 

The means and standard deviations of albedos and cloud 
top temperatures of each 10 x 10 pixel array for the weather 
types listed in Tables VIII and IX were calculated. Means 
and standard deviations of albedos and cloud top tempera- 
tures were also calculated for the 8x8, 6x6, and 4x4 

pixel arrays centered over the surface station. The 8x8, 
6x6, and 4x4 pixel arrays are equal to 36 nmi x 36 nmi, 
27 nmi x 27 nmi, and 22 nmi x 22 nmi at 30°N respectively. 
Variation of the digital satellite areal coverage is used to 
investigate the differences in the statistics due to the 
chosen resolution size. 

The data sets in Tables VIII and IX are represented, 
first, by the mean and standard deviation of the resolution 
cell means and standard deviations. Second, these data sets 



49 



are represented by the distributions of the mean cloud top 
temperatures and albedos where the mean cloud top tempera- 
tures are sorted into ten Kelvin (K) intervals and the mean 
albedos are sorted into 0.10 intervals. These representa- 
tive statistics and distributions are calculated for the 
four pixel array sizes. 

The statistical and distribution results for differing 
resolution sizes, bi-spectral threshold specification of 
precipitation, and separation of light from moderate/heavy 
precipitation are discussed in Chapter IV. 
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I V . RESDLT S 

A. INTRODUCTION 

The figures presented in this chapter display the dis- 
tributions of the grand means of the resolution cell means 
of albedos and cloud top temperatures for the data sets 
listed in Tables VIII and IX for the four array sizes. The 
mean values are sorted into ten Kelvin intervals and 0.10 
estimated albedo intervals. 

B. RESOLUTION 

The effect of satellite resolution in representing gen- 
eral precipitation (category 1) and no-precipitation over- 
cast (category 2A) data are explored for four resolution 
sizes. The four sizes are 10 x 10, 8 x 8, 6 x 6, and 4x4 

and are approximately equal to areas of 2025 nmi^, 1296 
nmi2 , 729 nmi^, and 484 nii^ at 30<*N respectively. 

The general precipitation (category 1) and no-precipita- 
tion overcast (category 2A) data were chosen for study 
because, while they represent two different weather condi- 
tions, their albedo and cloud top temperature distributions 
have the largest amount cf overlap when compared to any 
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other pair of precipitation versus no-precipitation data 
sets. The possibility arises that statistical differences 
in the four resolution sizes might be sufficient or comple- 
ment other information in delineating these two weather 
conditions. 

1 • Prec ipitation Data 

The general precipitation (category 1) data are com- 
prised of 329 overcast ceiling reports (61%), 205 broken 

ceiling reports (38%), and 4 scattered ceiling reports (1%). 

a. Mean Statistics 

The precipitation data (category 1) differences 
between the means of the cell means visual and infrared 10 x 
10 and 4 X 4 array sizes are 0.035 and 2.2K, respectively 
(Table X) . The trend of the mean of the means is toward 
higher albedo values and colder cloud top temperatures with 
the decreasing area or array size. The standard deviations 
of the means similiarly show an increase in the albedo, 
0.016, and cloud top temperature, 1.0K, • from the 10 x 10 

array size to the 4x4 array size. 

b. Standard Deviation Statistics 

The standard deviation statistics display the 
opposite trend with decreasing area as the mean statistics. 



52 



TABLE X 

Precipitation Data Statistics for Four Array Sizes 
Overcast, Broken, and Scattered Ceilings 







12. i 10. 


8x8 


6x6 


4x4 


Mean of 
Means 


(VIS) 


. 57 9 


. 591 


.603 


.614 




(IR) 


253. OK 


252. 2K 


251 . 4K 


250. 8K 




(-20OC) 


(-2 1°C) 


(-22°C) 


(-220C) 


Standard 

Deviations 


(VIS) 


. 21 1 


. 214 


.219 


. 227 


of Means 


(IR) 


20, IK 


2 1.0K 


21 .3K 


21. 7K 


Mean of 
St an dard 


(VIS) 


. 173 


. 161 


, 144 


. 1 24 


De viations 


(IR) 


1 0.4K 


9. IK 


7.6K 


5 .6K 


St an dard 
De viation 
of t he 
Standard 


(VIS) 


. 078 


. 078 


.075 


.072 


Deviations (IR) 


7. OK 


6.7K 


6 .OK 


4. 9K 



The means and standard deviations of 


the 


standard 


deviations 


decrease in the 


visual and infrared 


values 


with 


decreasing 


ar ea (Table X) . 


The differences between 


the 


10 X 


10 and 4 X 


4 array sizes 


visual and infrared 


means 


of the 


standard 



deviations are 0.049 and 4.8K, respectively, and the stan- 
dard deviations of the standard deviations are 0.006 and 
2. IK, respectively. 
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c. Distribution Discussion 

Tbs distributions of the precipitation data are 
shown in Figs. 10, 11, 12, and 13. There is a discernible 

shift toward higher albedcs and colder cloud top tempera- 
tures of the 2% and 3% frequency isopleth with decreasing 
array size. This upward shift is also reflected in the mean 
of the means (Table X). The appearance of the 5 % frequency 
isopleth in the 6x6 and 4 x 4 array sizes at high albedos 
and cold cloud tcp temperatures highlights the shift. 




Figure 10. Precipitation Data for 10 x 10 Array Size (The 
mean, .579 and 253. OK, interval is boxed. The 
2 % and 3 % frequencies are for 11 and 16 
occurrences, respectively.) 
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Figure 11. Precipitation Data for 8x8 Array Size (The 

mean, .591 and 252. 2K, interval is boxed. The 
2 % and 3% frequencies are for 11 and 16 
occurrences, respectively.) 




Figure 12. Precipitation Data for 6x6 Array Size (The 

mean. .603 and 251. 4K, interval is boxed. The 
2 %, 3X, and 5% frequencies are for 11, 16, and 
27 occurrences, respectively.) 
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Figure 13. Precipitation Data for 4x4 Array size (The 

mean, .614 and 250. 8K, interval is boxed. The 
2 %, 3%, and 5 % frequencies are for 11, 16, and 
27 occurrences, respectively.) 



The distributions of the 10 x 10 and 8x8 array 
sizes (Figs. 10 and 11) are unimodal while the 6x6 and 4 x 
4 array sizes (Figs. 12 and 13) appear to be more bimodal. 
The four array sizes were tested for a Gaussian distribution 
with the Chi-square test and all failed ar any confidence 
level. ■ Therefore, differences in the four resolutions can- 
not be adequately tested by well defined statistical methods 
based on an assumed normal distribution. 

The similiarities between the 10 x 10 and 8x8 
array sizes (Figs. 10 and 11) and the 6x6 and 4x4 array 



56 



sizes (Figs. 12 and 13) are further illustrated in Fig. 14. 
A diagonal cut is plotted for each of the four array sizes 
where the lines plotted are shown as a dashed boxed area in 
Figs. 10-13. The diagonal cut reveals the close agreement 
between the 6x6 and 4x4 array sizes along the line. The 
10 X 10 and 8x8 array' size lines follow the same general 
trend but do not coincide as closely as the 6x6 and 4x4 
array size lines. 

The chosen diagonal line results in the 8 x 8 
array size distribution appearing more smoothed than the 10 
X 10 (Fig. 14), as there is no relative minima at the 4.0 
interval for the 8x8 array size. The 10 x 10 array size, 
with the greater areal extent and therefore more averaging 
of differing clouds and clear areas, is expected to possess 
the •'smoothest” appearance, the lowest number of relative 
maxima and mirima of the four array sizes. However, the 8 x 
8 array size actually displays the fewest relative maxima 
and minima along the chosen diagonal line. The smoother 8 x 
8 array size cannot be explained in terms of significant 
differences in the number of cases of different ceiling 
types or different precipitation types occurring in interval 
4.0 between the four array sizes. Quite simply, the 
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smoother 8x8 array size apparently results from the sort 
ing intervals chosen for the distributions. 




CLOUD TOP TEMPERATURE/ ALBEDO INTERVALS 



Figure 14. Precipitation Array Size Distributions Along 

Diagonal Line (The line rep^resents the 10 x 10, 
dotted line the 8x8, dashed line the 6x6, 
and dash dot line the 4x4 array size.) 



The elongated shape of all four precipitation 



distribution array sizes (Figs. 10, 11, 12, and 13) reveal 



the variaticn in the areal amount of cloudiness and precipi- 



tation. The distributions range from high albedos and cold 



cloud top temperatures (indicative of satellite fields of 
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view filled with precipitating clouds) to low albedos and 
warm cloud top temperatures (indicative of satellite fields 
of view partially filled with precipitating clouds) . The 
visual and infrared satellite data distributions in Figs. 
10, 11, 12, and 13 agree with the elongated shapes of Platt 
(1981) for his cloud classifications of cumulus, frontal, 
and jetstream cirrus clcuds and agree with Coakley and 
Bretherton (1S82) for their general clouds present in a 1000 
lcm2 Pacific Ocean area, 
d. Summary 

A satellite field of view filled with a precipi- 
tating cloud is expected to have higher albedo and colder 
cloud top temperature values than a partially filled field 
of view. Additionally, the filled field of view would have 
a more uniform texture, as reflected in variance or standard 
deviation values, than a partially filled field of view. 
The statistics discussed in this study confirm these expec- 
tations for this data set. As the array size or field of 
view is decreased, the mean statistics increase while the 
standard deviation statistics decrease (Table X) . 

Fcr the precipitation data (category 1) , there 
are significant differences between the four resolution 
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sizes. These differences are reflected in the upward trend 
in albedos and cclder cloud top temperatures of the mean and 
the standard deviation of the means with decreasing area or 
array size (Table X) . The reverse trend is found in the 
mean and standard deviation of the standard deviations. 

The distributions have significant differences 
also. The relatively coarse resolution 10 x 10 and 8x8 
array sizes have a unimodal distribution while the rela- 
tively fine resolution 6 x 6 and 4 x 4 have a bimodal 
distribution. 

The four array sizes discussed vary in their 
statistics and distributions in representing the precipita- 
tion (category 1) data. The choice of satellite resolution 
for representation of the precipitation data will influence 
comparison of these data with other data. Therefore, for 
the remainder of this study the precipitation data for all 
classifications will be discussed using both the 10 x 10 and 
4x4 array size. , 

One additional topic to explore is that a number 
of cases with albedos less than 0.40 appear in the precipi- 
tation data (category 1) in all four size distributions 

(Figs. 10, 11, 12, and 13). These low albedo values suggest 
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the possibility that there might be a consistent low bias in 
the normalization scheme. But a closer look at the 
individual reports with albedos less than 0.40 show no 
pattern involved with either the GMT hour or the longitude 
or latitude of these station reports. Further analysis of 
these low albedo precipitation reports are discussed in the 
light precipitation section (IV. E.). 

2* No -pr e cip itatio n Overcast D ata 
a. Mean Statistics 

For the no-precipitation overcast cases (cat- 
egory 2A ) f differences between the means of the means visual 
and infrared 10 x 10 and 4x4 array sizes are 0.011 and 0.3 
K, respectively (Table XI). The 0.3K infrared difference is 
within the 0.5K noise level of the VIS3R infrared sensor. 
The trend of the visual mean of the means is upward with the 
decreasing array size. The standard deviations of the means 
show an increase in the albedo, 0.014, and cloud top temper- 
ature, 0.8K, from the 10 x 10 to 4 x 4 array sizes. Once 

again, both mean statistics have an upward trend with 

decreasing array size. 
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TABLE XI 

Nc-Precipitation Data Statistics for Four Array Sizes 





Overcast 


Ceiling 








11 i 11 


8x8 


6x6 


1x4 


Mean of (VIS) 


. 407 


. 411 


.415 


.418 


Means 

(IR) 


272. 6K 


272. 5K 


272. 4K 


272. 3K 




(-1°C) 


(- 1«C) 


(-10C) 


(-10C) 


Standard (VIS) 


. 205 


. 208 


.213 


. 219 


Deviations 










of Means (IR) 


1 7. 8K 


1 8. IK 


18 . 3K 


18. 6K 



Mean of (VIS) 


. 122 


. 113 


.102 


. 087 


St an dard 
De viations (IR) 


5.4K 


4.7K 


3.9K 


2.9K 


St an dard 


Deviation (VIS) 


. 056 


. 054 


.052 


.049 


of the 
St an dard 
Deviations (IR) 


4. 8K 


4.4K 


3 . 9K 


3. IK 



b. Standard Deviation Statistics 

Conversely, the standard deviation statistics 
have a downward trend with decreasing area size. The dif- 
ferences between the 10 x 10 and 4x4 array visual and 
infrared means of the standard deviations are 0.035 and 
2.5K, respectively. The differences in the standard devia- 
tions of the standard deviations are 0.007 in the visual and 
1.7K in the infrared values. 
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c. Distribution Discussion 

The distributions of the no-precipitation over- 
cast data are shown in Figs, 15, 16, 17, and 18. These dis- 

tributions are quite different from the precipitation 
distributions (Figs. 10, 11, 12, and 13). As expected, the 
no-precipitaticn overcast data are clustered at the low 
albedo and warm cloud top temperature values. The 5% fre- 
quency isopleths in Figs. 15, 16, and 17 show a grouping of 
the data at albedos ranging from 0.30 to 0.50 and cloud xop 
temperatures from 280K tc 290K. A bimodal distribution 
appears in the finer resolution 8x8, 6x6, and 4x4 

array sizes (Figs. 16, 17, and 18). This shifting of the 
no-precipitaticn overcast data into two relative maxima for 
the three smallest array sizes is the sols significant dif- 
ference in the four distributions. 

d. Summary 

The no-precipitation overcasr data display an 
upward trend in the two mean statistics with decreasing 
array size, although the 0, 3K mean of the means infrared 
difference is not significant. Conversely the two standard 
deviation statistics decrease with decreasing array size. 
These trends are consistent with the expected statistical 
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Figure 15. Nc-Precipitation Overcast Data for 10 x 10 Array 
Size (The mean, .407 and 212. SK., interval is 
boxed. The 2%, 3%, 5% and 1% frequencies are 
for 40, 59. 99, and 138 occurrences, 
respectively. ) 




CLOUD TOP TEMPERATURE (K) 

Figure 16. No-Precipitation Overcast Data for 8x8 Array 
Size (The mean, .411 and 272. 5K, interval is 
boxed. The 2%, 3%, 5% and 7% frequencies are 
for 40. 59, 99, and 138 occurrences, 
respectively. ) 
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CLOUD TOP TEMPERATURE (K) 

17. No-Precipitation Overcast Da- a for 6x6 Array 
Size (The mean, .415 and 272. 4K, interval is 
boxed. The 2%, 3%, and 5% frequencies are for 
40, 59, and 99 occurrences, respectively.) 




18. No-Precipitation Overcast Data for 4x4 Array 
Size (The mean, .418 and 272. 3K, interval is 
boxed. The 2^, 3%, and 5% frequencies are for 
40, 59, and 99 occurrences, respectively.) 
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trends discussed in the precipitation section (IV. B. I.)* 
However, in these data, there is a greater similiarity in 
the statistics for the four sizes because the ceilings are 
all overcast reports. 

The distributions (Figs. 15, 16, 17, and 18) 

provide visual confirmat icn of the similiarities between 
each of the four array sizes. with the exception of the 
second relative maxima at albedos of 0.30 zo 0.40 and cloud 
top temperatures of 260K tc 270K appearing in the 8x8, 6 x 
6, and 4x4 array sizes (Figs. 16, 17, and 18) , the four 
distributions are nearly identical. Because there are sta- 
tistical differences and distributional differences in the 
four array sizes, the 10 x 10 and 4x4 arrays sizes will be 
used to represent the two no-precipitation data categories. 

3* Precip itati on and Ho- pr e cipita tion compa riso n 

The possibility arises that statistical differences 
in the four resolution sizes might be sufficient, or at 
least complement other information, in delineating the pre- 
cipitation frcm the nc- precipita tion weather condition. The 
question then becomes, is there a statistic associated with 
variation of the array size within Tables X and XI which 
differentiates the precipitation reports from the no-precip- 
itation overcast reports? 
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If the satellite data processor can vary the 
resolution size, as was dene in this study by simply averag- 
ing different pixel array sizes, a trend in the visual and 
infrared data might be used to differentiate these two 
weather conditions. The most significant trend difference 

in the precipitation data (Table X) and no-precipitation 
overcast data (Table XI) occurs in the mean of the means. 
Recall that the mean of the means precipitation visual dif- 
ference between the 10 x 10 and 4x4 array sizes was 0.035 
while the no-precipitation overcast visual difference was 
0.011. Similiarly the infrared differences were 2.2K and 
0.3K for the precipitation and no-precipitation overcast, 
respectively. The precipitation data show a greater upward 
trend toward higher albedos and colder cloud top tempera- 
tures than the no-precipitation overcast with the finer sat- 
ellite resolution. 

Differentiation between these two data sets based on 
a comparison of the trend in the mean of rhe means are sug- 
gested by Tables X and XI. It must be emphasized that these 
tables are based on many reports and therefore reflect the 
most rypical values. Individual reports within a given 
interval should he studied to provide conclusive evidence as 
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to whether these statistics can be used on a few reports to 
differentiate precipitat icn from no-precipitation overcast 
re ports. 

C. PRECIPITATICN SPECIFICATION 

An essential difference between this specification study 
and most of those in the literature is that the distribu- 
tions of the precipitation and no-precipitation data sets 
are examined in detail to extract information about the 
probability of correct classifications. Only Lovejoy and 
Austin (1979) present their data distributions. The bi- 
spectral and life history method thresholds (Tables I and 
VII) refer to the typical or most common threshold values 
for precipitation, which is assumed to be eguivalenr to the 
mean of the means in this study. Therefore, the mean of the 
means can be coirpared to the threshold values in Tables I 
and VII. Additionally, a bi-spectral threshold can be pro- 
posed based on the distributions and with these distribu- 
tions the amount of overlap, or the percentage of correctly 
classified precipitation or no-precipitation cases, can be 
calculated. 
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and 



"I* 2Z6£cast Ceilings 

Are the precipitation overcast (category 1A) 
no-precipitation overcast (category 2A) data sets suffi- 
ciently separated to allow differentiation of the two popu- 
lations? If so r how much overlap is there between the two 
data sets? 

a. Mean of the Means 

The mean of the means statistics for the precip- 
itation overcast versus no -pr ecipitation overcast data. 
Table XII, show there is a .242 and 24. 4K difference between 
the two populations for the 10 x 10 array size and a .254 
and 25. 9K difference for the 4x4 array size. The respec- 

tive differences are greater than one standard deviation of 
the means of either of the two populations. 

If the two precipitation overcast array sizes 
are compared to the no- preci pitation overcast array sizes, 
the stronger trend in the mean of the means (Table XII) is 
seen in the precipitation overcast data. While the precipi- 
tation visual and infrared values vary by 0.023 and 1.8K 
between the 10 x 10 and 4x4 array sizes, the no- precipita- 
tion visual and infrared values vary by 0.011 and 0.3K. 
Both the mean cf the means and their trends can be used to 
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TABLE XII 

Precipitation Specification Overcast Ceilings 







Precipi- 

tation 


No -Preci- 
pitat ion 


Precipi- 

tation 


No-Preci 

pitation 






10 2 10 


10 X J[0 


it 2 i 


4x4 


Mean of 
Means 


(VIS) 

(IR) 


. 64 9 

248. 2K 
(-250C) 


. 407 

27 2. 6K 
(- IOC) 


.672 

246 . 4K 
(-27®C) 


.418 

272. 3K 
(-10C) 


Stan dard 
Deviations 
of Means 


(VIS) 

(IR) 


. 186 
1 8. 9K 


. 205 
1 7. 8K 


.201 
19. 8K 


.219 
18. 6K 



Mean of (VIS) 


. 152 


. 122 


.108 


. 087 


St an dard 
De viations (IR) 


8.9K 


5.4K 


4 .4K 


2.9K 


Standard 


Deviation (VIS) 


. 078 


. 056 


.069 


. 049 


of t he 

Standard 

De viations (IR) 


6. 3K 


4. 8K 


4. OK 


3 . IK 



differentiate precipitation from no-pr ecipiration for a 
large number of reports in a region similiar to this summer- 
time convective shower dominated area. 

b. Mean of the Standard Deviations 

The visual and infrared means of the standard 
deviations vary by 0.030 and 3.5K for rha 10 x 10 array size 
and by 0.021 and 1.5K for the 4 x 4 array size. In the 
infrared values, the no-precipitation mean of the standard 
deviations have a magnitude 60% of the precipitation values 
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and produce relatively large differences. In the visual 
values, the no-precip itat icn mean of the standard deviations 
have a magnitude 80% of the precipitation values. The rela- 
tively large differences (cnly the mean of the means have a 
larger difference) suggest the use of this statistic to dif- 
ferentiate precipitation overcast from no- precipitation 
overcast. 

c. standard Deviation of the Means 

The differences in the standard deviations of 
the means (Table XII) are nearly equal when comparing the 
two 10 X 10 array sizes and the two 4x4 array sizes. The 

visual differences are 0.019 for the 10 x 10 and 0.018 for 
the 4x4 array size. Similiarly, the infrared differences 
are 1.1K for the 10 x 10 and 1.2K for the 4x4 array size. 
The differences in the 4x4 array size (0.018 and 1.2K) are 
comparable to the relatively significant differences in the 
means of the standard deviations (.021 and 1.5K). However, 
the differences are not comparable in the 10 x 10 array 
size . 

d. Standard Deviation of the Standard Deviations 
The differences in the srandard deviations of 

the standard deviations (Table XII) for the two array sizes 
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are approximately equal also. The visual differences are 
0.022 for the 10 x 10 and 0.020 for the 4 x 4 array size. 
The infrared differences are 1.5K for the 10 x 10 and 0.9K 
for the 4x4 array size. Once again the differences in the 
4x4 array size (0.020 and 0.9K) are approximately equal to 
the differences in the means of the standard deviations 
(0.021 and 1.5K), particularly in the visual value, 
e. Distribution Discussion 

The distributions for the 10 x 10 array size 
precipitation overcast (Fig, 19) and no-precipitation over- 
cast (Fig. 15) and the 4x4 array size precipitation over- 
cast (Fig. 20) and no- precipi tation overcasr ((Fig. 18) 
allow visual confirmation of the degree of separation of 
these two populations. These figures verify the separation 
between the the occurrence maxima of the two populations 
while showing that there is overlap of some of the values in 
the two populations. 

The appearance of a bimodal distribution in the 
relatively fine resolution 4 x 4 array size precipitation 
overcast (category 1A) and no- precipitation overcast (cat- 
egory 2A) in Figs. 20 and 18 cannot be explained in terms of 
the precipitation categories (Table IX) or ceiling cover. 
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Figure 20. Precipitation Overcast Data for 4x4 Array Size 
(The mean, .672 and 246. 4K, interval is boxed. 
The 2%, 3%, 555 and 7% frequencies are for 7, 10, 
16, and 23 occurrences, respectively.) 
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An alternate explanation might be that these two maxima 
reflect different synoptic signatures. Four weak frontal 
systems impact this data set region during August 1979 and 
cause a change in the cloudiness and precipitation pattern 
which is normally produced by daytime heating. Recall that 
the Lovejoy and Austin (1979) data set also showed bimodal 
distributions for the cumulus no-rain reports. The bimodal 

distributions in this study may not be due to a synoptic 
signature. Nonetheless, this possibility should be 
investigated. 

f. Precipitation Probabilities 

Precipitation probabilities (Figs. 21 and 22) 
were computed from the precipitation overcast (category 1A) 
and the no-precipitation overcast (category 2A) data for the 
10 X 10 (Figs. 19 and 15) and the 4x4 (Figs. 20 and 18) 
array sizes. Estimated albedos greater than 1.00 were not 
included in these probabilities as they accounted for only 
two and three no-precipitation overcast reports and four and 
five precipitation overcast reports in the 10 x 10 size and 
4x4 array size, respectively. The two probability figures 
indicate that the 5056 probability of precipitation is not a 
simple function of mean albedo and mean cloud top 
temperature. 
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Figure 21. Precipitation Overcast Data Probability 10 x 10 
Array Size 




Figure 22. Precipitation Overcast Data Probability 4x4 
Array Size ^ 



The 50% probability line (Fig. 21) shows the 
precipitation at low clcud top temperatures (270K-290K) 
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occurs at high albedos (0.80-1. 00) and at cold cloud top 
temperatures (210K-230K) occurs at relatively low albedos 
(0.40-0.60). One exception occurs in the 0.60-0.80 albedo 
and 210K-230K cloud top temperature interval and represents 
26 precipitaticn reports cf 64 total reports. The 100% pre- 
cipitation probability at 0.00-0.20 albedo and 230K-250K 
results from two precipitation reports. Two reports in a 
0.20 albedc and 20K cloud top temperature interval are not a 
sufficient number of reports to be a significant indication 
of a high probability of precipitation. The 50% precipita- 
tion probability line (Fig. 22) in the 4x4 array size data 
shows the same general trend of decreasing albedo with 
decreasing cloud top temperatures. In this fine resolution 
data (Fig. 22 ) , there is an upturning of the 50% probability 
line at the coldest cloud top temperatures, 210K-230K. 
There are 73 no-precipitaticn overcast and 80 precipitation 
overcast reports in the 4x4 array size dara in the 
210K-230K interval so the upturning is nor the result of 
lack of data. The appearance, once again, of a greater than 
50% probability cf precipitation at low albedos, 0.00-0.20, 
between 210K-250K results from a total of four reports (one 
precipitation report of two total reports in the 230K-250K 
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interval and two precipitation reports of two reports in the 
210K-230K) and is not a significant indication of high pre- 
cipitation probability. 

If a straight line is drawn to represent the 50% 
probability line, the lines for the two array sizes are 
nearly coincident until the 230 K cloud top temperature is 
reached. For the purposes of this study, a simple linear 
function bi-spectral precipitation threshold based on a 50% 
probability of precipitation can be approximately defined as 
extending from 1.00 albedo and 290K cloud top temperature to 
0.60 albedo and 210K cloud top temperature. 

Comparison of this linear bi-spectral threshold 
with the Muench and Keegan (1979) threshold (Fig. 2), shows 
the proposed thr^hold has an albedo approximately 0.10 
smaller at corresponding cloud top temperatures. The 50% 
probability asymptote at the warm cloud top temperatures 
shown in the Muench and Keegan (1979) results (Fig. 2), are 
not shown in Figs. 21 and 22 due to lack of reporrs in these 
values. The lower albedo values for this linear bi-spectral 
threshold may be associated with the dominance of convective 
precipitation (426 of 538 reports) in the precipitation 
data. The effect of convective precipitation reports on the 
satellite albedo values is discussed in Section IV. D. 
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g. Summary 

The two significant statistics for precipitation 
specification of overcast ceiling reports for both array 
sizes are the mean of the means and the mean of the standard 
deviations (Table XII). The probabilites in Figs. 21 and 22 
make use of the mean of the means only. Inclusion of the 
mean of the standard deviations in the precipitation and 
no-precipitaticn distribution plots may more distinctly 
define the two weather data types. The question is, how to 
graphically display four variables (i.e. four dimensions) in 
one plot? 

One solution is to find a three-dimensional plot 
that involves the four variables. The most straight forward 
approach is tc define the 50% probability for the mean of 
the means (Figs. 21 and 22) in terms of a surface. A plane 
would be the simplest surface choice. The equation of a 
line perpendicular to the 50% probability plane intersecting 
its midpoint in the distribution planes shown would then be 
calculated. All of the precipitation and no-precipitation 
points would be projected onto the line and the line would 
become the x-axis in a new plot. Thus this x-axis reflects 
the visual and infrared mean of the means. The visual and 
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infrared mean of the standard deviations would define the 



y-axis and z-axis, respectively. 

The defined plot for the precipitation and no- 
precipitation data would show the relative dependence of the 
data on the four variables. If the plots produced distinct 
groups for the precipitation and no-precipitation data, new 
precipitation probabilities would be calculated. The varia- 
tion with satellite resolution size of the precipitation 
probabilities would then have to be reconsidered with the 
new data. The data processing described is beyond the scops 
of this particular research effort and is reco mm mended for 
future investigation. 

Precipitation specification for overcast ceil- 
ings can he delineated by the values for the mean of the 
means and the mean of the standard deviations for any array 
size discussed in this study. The mean of the mean values 
when used with Figs. 21 and 22 will indicate the probability 
of precipitation, given a similiar time of year and climato- 
logical area. A simple linear bi-spectral threshold, based 
on a 50% probability of precipitation, is defined approxi- 
mately as extending from a 1.00 albedo and a 290K cloud top 
temperature to a 0.60 albedo and a 210K cloud top 
temperature. 
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are the precipita- 



2* Overcast an d Broken C ei l inq s 

Once again the question arises, 
tion overcast and broken (category IB) and no-precipitation 
overcast and broken (category 2B) data sets sufficiently 
separated to allow differentiation of the two populations? 
If sc, how much overlap is there between the two data sets? 
a. Mean of the Means 

The mean of the means statistics for the precip- 
itation overcast and broken (category 1B) versus no-precipi- 
tation overcast and broken (category 2B) data. Table XIII, 
show there is a .309 and 27. OK difference between the two 
populations for the 10 x 10 pixel size and a .339 and 29. IK 
difference for the 4x4 pixel size. The respective differ- 
ences are approximately equal to one and one-half standard 
deviations of the means of either of the two populations. 

For these data, the differences in the trends of 
the mean of the means between the 10 x 10 and 4 x 4 array 
sizes for the precLpitaticn overcast and broken reports are 
more dramatic than for the precipitation overcast reports. 
Thera is a 0.036 visual and a 2.2K infrared difference in 
the precipitation overcast and broken data and only a 0.006 
visual and a 0.1K in the no-precipitation overcast and 
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TABLE XIII 

Precipitation Specification Overcast and Broken Ceilings 



Mean of (VIS) 
Means 

(IR) 



Standard (VIS) 
De viations 
of Means (IR) 



Mean of (VIS) 
Standard 
Deviations (IR) 



Standard 
Deviation (VIS) 
of the 
St an dard 
Deviations (IR) 



broken data between 



Therefore, variation 



Preci pi- 
tation 


No -Preci 
pitation 


10 X 10 


10 X 10 


. 580 


. 271 


252. 8K 
-200 C) 


27 9.8K 
(+70C) 


. 210 


. 189 


20.5K 


1 6.6K 


. 173 


. 128 


1 0.4K 


5.9K 


. 078 


. 064 


7. OK 


5.7k 



10 and 
re solution 
be used. 



Precipi- 

tation 


No-Preci- 

pitation 


4x4 


i£ X i 


.616 


. 277 


250. 6K 
(-230C) 


279. 7K 
( + 70 C) 


.225 


. 207 


21. 5K 


17 . 8K 


.124 


. 094 


5. 6K 


3. 2K 


.072 


. 057 


4 .9K 


3. 7K 


X 4 arr 


ay sizes. 


size and 


the trend 



as well as the mean 



the 10 X 
of the 

for the mean of the means can 



of the means value itself, in delineating precipitation from 
no-precipitation for these overcast and broken ceiling 
reports. 

b. Mean of the Standard Deviations 

The values of the mean of the standard devia- 
tions are important. The 0.045 albedo and 4.5K cloud top 
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temperature difference between the two data sets for the 10 
X 10 array size are significant as well as the 0.030 and 
2.4K difference for the 4x4 array size. 

c. standard Deviations of the Means 

The differences in the standard deviations of 
the means are nearly equal when comparing the two 10 x 10 
array sizes and the two 4x4 array sizes. The visual dif- 
ferences are 0.021 for the 10 x 10 and 0.018 for the 4x4 
array size. Similiarly, the infrared differences are 3.9K 
for the 10 x 10 and 3.7K for the 4x4 array size. The vis- 
ual differences are 40% less than the statistically signifi- 
cant visual mean of the standard deviation differences and 
the infrared differences are comparable to the infrared mean 
of the standard deviation differences. 

d. Standard Deviation of the Standard Deviations 
The differences in the standard deviations of 

the standard deviations are approximately equal. The visual 
differences are 0.014 for the 10 x 10 and 0.015 for the 4 x 
4 array size. The infrared differences are 1.3K for the 10 
X 10 and 1.2K for the 4x4 array size. Both the visual and 
infrared differences are 5055 less than the respective sig- 
nificant differences in the mean of the standard deviations. 
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e. Distribution Discussion 

The distributions for the 10 x 10 array size 
precipitation overcast and broken (Fig. 23) and no-precipi- 
taticn overcast and broken (Fig. 24) and the 4 x 4 array 
size precipitation overcast and broken (Fig. 25) and no-pre- 
cipitation overcast and broken (Fig. 26) allow visual con- 
firmation of the relatively larger degree of separation 
between these two populations compared to the overcast ceil- 
ing data. These figures once again also show overlap 
between the two populations. 




Figure 23. Precipitation Overcast and Broken Data for 10 x 
10 Array (The mean, .580 and 252. 8K, interval is 
boxed. The 2J[ and 3% frequencies are for 11 and 
16 occurrences, respectively.) 
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Figure 24. Nc-Precipitat icn Overcast and Broken for 10 x 10 
Array (Tne mean, ,271 and 279. 8K, interval is 
boxed. The 2%., 2%, S%, and 7% frequencies are 
for 147, 221, 368, and 515 occurrences, 
respectively. ) 




Figure 25. Precipitation Overcast and Broken Data for 4x4 
Array (The mean, .616 and 250. 6K, interval is 
boxed. The 2%, 3%, and 5% frequencies are for 
11, 16, and 27 occurrences, respectively.) 
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Figure 26. No-Precipita t icn Overcast and Broken for 4x4 
Array (The mean, .277 and 279. 7K, interval is 
boxed. The 2 %, 355. 5%, 755, and 12% frequencies 
are for 147, 221, 368, 515, and 883 occurrences, 
respectively. ) 



f. Precipitation Probabilities 

Precipitation probabilities (Figs. 27 and 28) 
were computed from the precipitation overcast and broken 
(category IB) and the no-precipitation overcast and broken 
(category 2B) cata for the 10 x 10 (Figs. 23 and 24) and the 
4x4 (Figs. 25 and 26) array sizes. As in the overcast 

ceiling reports, there is a greater than 50% probability of 
precipitation at the low albedo values, 0.00-0.20. In Fig. 
27, the 75% probability between 230K-250K interval results 
from three precipitation reports of four total reports. The 



I 
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overcast ceiling reports account for two of the 
precipitation reports. Therefore, one precipitation and one 
no-precipitaticn broken report have been added to the inter- 
val. Similiarly in Fig. 28, the 67% probability in the 
210K-230K interval results from two precipitation reports of 
three total reports. The overcast ceiling reports account 
for two of the precipitation reports. Thus, one no-precipi- 
tation report has been added to the interval. A few reports 
in these low albedo and cold cloud top temperature intervals 
are producing misleadingly high precipitation probabilities. 
The 50% probability line is nearly constant, at a 0.80 
albedo, at all cloud top temperatures for the overcast and 
broken ceiling data. 

The same treatment for displaying the four vari- 
ables in a three dimensional diagram suggested in the previ- 
ous section (IV.C. 1.), is recommended for these data as 
well. The new displays then could be used to calculate pre- 
cipitation probabilities that might allow a more accurate 
bi-spectral threshold specification of precipitation, 
g. Summmary 

As for the overcast ceilings, precipitation 
specification of the overcast and broken ceilings can be 
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Figure 28. Precipitation Data Probability 4x4 Array Size 
(Overcast and Eroksn Ceilings) 



delineated by the values for the mean of the means and the 
mean of the standard deviations for any array size discussed 
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in this study. The mean of the means values when used with 
Figs. 27 and 28 will indicate the probability of precipita- 
tion, given a similiar time of year and climatological area. 
A very simple threshold fcr these data is actually dependent 
upon the visual data value being greater than 0.80. 

D. CONVECTIVE VERSUS CONTINUOUS PRECIPITATION 

The physical processes involved in convective precipita- 
tion, or showers, are different from those processes usually 
involved in continuous precipitation. The inherent differ- 
ences in the processes might result in a statistical separa- 
tion in the thresholds between the two types of 
precipitation. For this data set, 426 of the 538 cases are 
classified as convective by the surface weather report, 
with the data set consisting of southeastern United States 
stations in August, the dominance of the precipitation cases 
by convective reports is not surprising. 

Wean Stat ist ics 

The statistics for the convective and continuous 
precipitation (Table XIV) indicate differences in the means 
of the cell means of .013 and 5.7K for the 10 x 10 array 
size and differences of .015 and 6.8K for the 4x4 array 
size. As mentioned in the life history methods, convective 
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precipitation is expected to be associated with relatively 
higher albedos and colder cloud top temperatures. For these 
data, the convective precipitation are associated with 



colder cloud top 


temperatures , 


resulting 


from the 


cumulonim- 


bus clouds. 


TABLE 


XIV 






Continuous versus Convective 


Precipitation Specification 




Contin- 


Co nvec- 


Contin- 


Convec- 




uous 


t ive 


uous 


tive 




10 X 10 


10 X 10 


4x4 


4x4 


Mean of (VIS) 

Means 


.58 9 


. 576 


.602 


. 6 17 


(IR) 


257.5K 


25 1.8K 


256. 2K 


249. 4K 


(-160C) ( 


-2 1°C) 


(- 170C) 


(-24«C) 


Standard (VIS) 
De viations 


. 201 


. 214 • 


.221 


. 228 


of Means (IR) 


1 9. 4K 


2 0.8K 


20 . IK 


21 .9K 


Mean of (VIS) 


. 147 


. 180 


.099 


.131 


Standard 
Deviations (IR) 


7. 6K 


1 1. IK 


3.8K 


6. OK 


Standard 










Deviation (VIS) 
of the 


.075 


. 078 


.058 


.075 


Standard 
Deviations (IR) 


5. 3K 


7.2K 


3. 2K 


5. 2K 



The statistics for the convective and continuous 
precipitation (Table XIV) suggest that for visual satellite 
values, “he mean of the standard deviations is the best sta- 
tistic for differentiation of these two data types, while 
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for infrared values. 



the mean of the means is the best 



statistic. The use of two different statistical types to 
qualitatively specify continuous precipitation from convec- 
tive precipitation is unique to this pair of convective ver- 
sus continuous precipitation data comparison. 

The infrared mean of the cell means show that the 
convective precipitation has colder cloud top temperatures 
due to the greater vertical development of the clouds. The 
greater vertical variation of the convective precipitation 
clouds is seen in the visual "texture" or visual mean of the 
standard deviation statistics. This "texture" difference 
also appears in the infrared values, but it is not as sig- 
nificant as the differences in the infrared mean of the 
means. The lack of a large difference in the visual mean of 
the means between the convective and continuous precipita- 
tion reflects the averaging of open areas and cumulus clouds 
associated with convective precipitation. Note that for the 
fine resolution 4x4 array size, the visual mean of the 
means is larger for the convective than the continuous 
precipitation. 

The visual means of the standard deviations have a 
difference of 0.033 for the 10 x 10 and 0.032 for the 4x4 
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array size. Contrast these difference to the visual mean of 
the mean differences of 0.013 and 0.015 for the 10 x 10 and 
4x4 array size, respectively. The infrared mean of the 
means have a difference of 5.7 k for the 10 x 10 and 6.8K for 
the 4x4 array size. These differences are larger than the 
infrared mean of the standard deviations differences of 3.5K 
for the 10 x 10 and 2.2K for the 4x4 array size. 

2* St anda rd De v iat ion Stati sties 

The standard deviation of the means and the standard 
deviation of the standard deviations display quite similiar 
and relatively insignificant differences. The visual stan- 
dard deviation of the means vary by 0.013 in the 10 x 10 and 
0.007 in the 4x4 with the infrared values varying by 1.4K 
for the 10 x 10 and 1.8K for the 4x4 array sizes. The 
visual and infrared 10 x 10 and 4x4 array size differences 
in the standard deviations of the standard deviations are 
0. 003, 0.017, 1.9K, and 2. OK, respectively. 

3* Distribu tio n Discu ssion 

The distributions of the convective and continuous 
precipitation for the 10 x 10 array size (Figs. 29 and 30) 
and for the 4x4 array size (Figs. 31 and 32) demonstrate 
the similiarity of the distributions. This similiarity 
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disallows quantitative separation of the convective and con- 



tinuous precipitation. 




Figure 29. Convective Pr-ecipita tion Data for 10 x 10 Array 
Size (The mean, , 576 and 251. 8K, interval is 
boxed. The 25?, 3%, and 5% frequencies are for 
9, 13, and 21 occurrences, respectively.) 



4 . Su mmary 



Defining the qualitative delineation of convective 



versus continuous precipitation in terms of the shift toward 



colder cloud top temperatures for the convective precipita 



tion cases is one available delimiter. The second qualita- 



tive separation is in terms of the shift toward higher 



albedo values of the mean of the standard deviations for the 
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Figure 30. Continuous Precipitation Data for 10 x 10 Array 
Size (The mean, .589 and 257. 5K, interval is 
• boxed. The 2%, 3%, and 5% frequencies are for 

2, 3, and 6 occurrences, respectively.) 




CLOUD TOP TEMPERATURE (K) 

Figure 31. Convective Precipitation Data for 4x4 Array 
Size (The mean, .617 and 249. 4 K, interval is 
boxed. The 25E, 3%, and 5% frequencies are for" 
9, 13, and 21 occurrences, respectively.) 
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CLOUD TOP TEMPERATURE (K) 

Figure 32. Ccntinuous Precipitation Data for 4x4 Array 
Size (The mean, .602 and 256. 2K, interval is 
boxed. The 2%, 3%, and 5% frequencies are for 
2, 3, and 6 occurrences, respectively.) 



convecrive precipitation cases as compared to the continous 
precipitation cases. 

The convective precipitation cloud top temperature 
mean of 251. 8K (-210C) ccmpares quite well with those asso- 
ciated with life history methods (Table VII) . The continu- 
ous precipitation cloud top temperature of 257. 5K (-16®C) 
compares well with the -16®C threshold of Wylie (1979) (a 
life history threshold for Montreal, Canada, in summertime) 
and with the -12®C bi-spectral threshold of Muench and Kee- 
gan (1979). These research efforts focused on mid-latitude 
precipitation associated with fronts. 
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E. INTENSITY 

As there is little statistical difference between the 

convective and continuous precipitation data, the gualita- 

\ 

tive analysis of precipitation can be based strictly on all 
light or all moderate/heavy precipitation. The light pre- 
cipitation cases are expected to have a lower mean albedo 
and warmer cloud top temperature when compared ro the moder- 
ate/heavy precipitation cases. 

^ • Stat istical Disc ussio n 

The statistics fcr the two precipitation intensi- 
ties, Table XV, show 0.100 and 8.6K diffarences for the 10 x 
10 array size and .117 and 8. IK differences for t-he 4 x 4 
array size in the means of the means. The differences are 
approximately egual to one-half a standard deviation, an 
insufficient separation for specification of these two pre- 
cipitation classifications. The remaining three sets of 
statistics (Table XV) show nearly identical values in the 
light and moderate/heavy precipitation for the two array 
si zes. 



95 



TABLE XV 

Light versus Moderate/Heavy Precipitation Specification 





Light 


Moderate/ 


Light 


Moderate/ 




He avy 


Heavy 




11 X 10 


10 X 10 


4x4 


4x4 


Mean of (VIS) 

Means 


. 56 5 


. 665 


.598 


.715 


(IR) 


254. 2K 


24 5.6K 


251 . 9K 


243 . 8K 


(- 19®C) 


(-2 8®C) ( 


-210C) 


(-29«C) 


Standard (VIS) 

Deviations 


. 210 


. 201 


.226 


. 204 


of Means (IR) 


20. 5K 


2 0.4K 


21 .6K 


21 . 3K 


Mean of (VIS) 

Standard 


. 172 


. 177 


. 124 


. 129 


De viation s (IR) 


10. 2K 


1 1.9K 


5.5K 


6 . IK 


Standard 
Deviation (VIS) 
of the 


. 078 


. 081 


.072 


. 079 


St andard 
De viations (IR) 


7. OK 


7. OK 


4.9K 


5.2K 


2. Distribution Discussion 






The distributions of 


the light 


and mode 


rate/heavy 


precipitation for 


the 10 X 10 


(Figs. 33 


and 34) and the 4 x 


4 array size (Figs. 35 and 36) allow visual confi 


rmation of 


the overlap between the two 


in tensities 


. Basi 


cally with 


only 74 cases of 


mod era t e/heavy precipita 


tion, the 


frequency 


contours (Figs. 


34 and 36) 


reflect a 


noisy di 


str ibution 



created by a sorting of just two or three reports into a 
given interval. Additional intensity data, specifically 
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moderate/heavy prscipitation cases, are needed before a 
threshold delineating qualitative intensities can be 
proposed . 




Figure 33. Light Precipitation Data for 10 x 10 Array Size 
(The mean, .565 and 254. 2K, interval is boxed. 
The 2%, 3 %, and 5% frequencies are for 9, 14, 
and 23 occurrences, respectively.) 



For the 10 x 10 and 4x4 array sizes moderate/heavy 



precipitation (Figs. 34 and 36) , the area encompassed by the 



7% frequency isopleth at the higher albedos and colder cloud 



top temperatures gives an indication of a possible 



separation between the two intensities. However, the bimo- 



dal nature of the 10 x 10 distribution (Fig. 34) and the 
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Figure 34. Moderate/Heavy Precipitation Data for 10* x 10 

Array Size (The mean, .665 and 245. 6K, interval 
is boxed. The 3%, 5%, and 1% frequencies are 
for 2, 4, and 5 occurrences, respectively.) 




CLOUD TOP TEMPERATURE (K) 



Figure 35. Liaht Precipitation Data for 4x4 Array Size 
(The mean, .598 and 251. 9K, interval is boxed. 
The 2% and 3% frequencies are for 9 and 14 
occurrences, respectively.) 
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Figure 36. Moderate/Heav y Precipitation Data for 4x4 

Array size (T&e mean, ,715 and 243, 8K, interval 
is boxed. The 3%, S%, and 7% frequencies are 
for 2, 4, and 5 occurrences, respectively.) 



three relative maxima of the 4x4 distribution (Fig. 36) 



suggest that there will be a significant amount of overlap 



between these two qualitative intensities. Because the sta- 



tistics for the two rainfall rates are nearly identical 



(Table XV) , with the exception of the mean of the means, the 



method discussed in the precipitation specification section 



(IV. c.) which involved the mean of the standard deviations 



to aid in further differentiation of two data types cannot 



be applied to these two intensities. The delineation will 



have to be derived from the mean of the mean statistics and 



their resultant distributions. 
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3. Summary 

There are slight manifestations of the expected 
downward shift in the albedos and the warmer cloud top temp- 
eratures for the light precipitation relative to the moder- 
ate/heavy precipitation. However, additional data are 
required before a threshold can be proposed to delineate 



light versus 


mode ra te/heav y 


preci pitation. 


The 


investigation 


of 


the 


addit ion al 


data should begin 


with 


cons ideraticn 


of 


the 


mean of 


the means and 


their 


di stribut ions. 












During 


the 


cell 


resolution discussion of Section 


IV.B.1., a significant 


number of 


low albedo values 


(less 


than 0.40 albedc) 


were 


noted. A 


comparison of the 


light 


precipitation 


data 


(Fig. 


33) and 


precipitation data 


(Fig. 



10) for the 10 x 10 array size reveals that 96% (105 of 109 
reports) are cl^sified as light precipitation. Further 
analysis shows that 79% (86 of 109 cases) are classified as 
light convective precipitation. (The decomposition of con- 
vective and continuous precipitation data into light and 
moderate/heavy are not shown.) These low albedo values are 
not unreasonable for light convective precipitation which 
can be produced by small, isolated cumulus congestus clouds. 
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V. SOM MARY AND CON C LOS 10 NS 
A. DATA PROCESSING SOMMAEY 

The data were first filtered for the correct time, 
1200-2000 GMT, to ensure reliable visual satellite values. 
The visual counts were converted to albedos according to the 
Muench and Keegan (1979) normalization scheme (see Appendix 
A) . This scheme corrects for Lambertian scattering as well 
as anisotropic cloud radiation. Sorting the data according 
to the Service-A reports (current weather and cloud group) 
provided seven precipitation groups: convective, continuous, 
light, moderate/heavy, general, overcast ceiling, and over- 
cast and broken ceiling precipitation and two no-precipita- 
tion groups: no-precipitation overcast and no-precipitation 

overcast and broken. In order to investigate the impact of 
satellite resolution on precipitation specification, the 
stated groups were subdivided into four sizes: 10 x 10 (45 x 

45 nmi) , 8 x 8 (36 x 36 n mi) , 6 x 6 (27 x 27 nmi) , and 4x4 
(22 X 22 nmi) pixel array sizes. The means of the albedos 
and cloud top temperatures for each array size for each 
weather or cloud group were calculated. The mean values 



101 



ware sorted into 10K cloud top temperature intervals and 
0.10 estimated albedo intervals to provide a visual descrip- 
tion of the distributions. The standard deviation of the 
means and the mean and standard deviation of the standard 
deviations for both albedc and cloud top temperature were 
also computed for each distribution. 

# 

B. STATISTICS SUMMARY 

The computed mean and standard deviation of the means 
and the mean and standard deviation of the standard devia- 
tions for both albedo and cloud top temperature were evalu- 
ated for specification of precipitation, specification of 
continuous versus convective precipitation, and specifica- 
tion of light versus moderate/heavy precipitation. 

For precipitation specification (Tables XII and XIII), 
the visual and infrared mean of the means produced the larg- 
est statistical differences between the precipitation and 
no-precipition data sets fcr both the overcast and the over- 
cast and broken ceiling reports. For the overcast ceiling 
reports (Table XII), the separation in the visual and infra- 
red means of the means were one standard deviation. For the 
overcast and broken ceiling reports (Table XIII), the sepa- 
ration in the visual and infrared means of the means were 
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one and one-half standard deviations. The mean of nhe stan- 
dard deviations and the trend of the mean of the means also 
yielded a large difference between the precipitation and 
no-precipitarion data sets. 

For specification of continuous versus convective pre- 
cipitation (Table XIV), the infrared mean of the means and 
the visual mean of the standard deviations yielded the larg- 
est statistical differences. The convective precipitation 
infrared mean of the means was 6K-7K colder than the contin- 
uous precipitation value and the convective visual mean of 
the standard deviations was 20% greater than the continuous 
precipitation value. These two statistics, however, were 
only indications of possible separation as the two data dis- 
tributions displayed significant overlap. 

For light versus moderate/heavy precipitation specifica- 
tion (Table XV), the mean of the means yielded the largest 
statistical differences. However, the relatively small 
quantitative differences (one-half standard deviation) in 
the visual and infrared values only provided an indication 
of the upward shift in albedo and the downward shift in 
cloud top temperature for the moderate/hea vy precipitation 
relative to the light precipitation. 
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C. PRECIPITATION PROBABILITY SUMMARY 

For the overcast ceiling reports, the 10 x 10 array size 
(Fig. 21) 50 % probability line decreased linearly to lower 

albedos with colder cloud top temperatures, while the 4x4 
array size (Fig. 22) 50% probability line decreased linearly 

to approximately 240K and then increased linearly at colder 
cloud top temperatures. For the overcast and broken ceiling 
reports, the 10 x 10 array size (Fig. 27) 50% probability 

line was constant at 0.80 albedo while the 4x4 array size 
(Fig. 28) 50% probability line gradually decreased linearly 
to approximately 240K and then gradually increased linearly 
at cclder cloud top temperatures. 

D. DATA DISTRIBUTION SUMMARY 

The distributions of all seven precipitation classes and 
both no-precipitation classes were non-Gaussian. This 
result disagrees with the assertion of Lovejoy and Austin 
(1979) that their cumulus rain, non-cumulus rain, and non- 
cumulus no-rain data distributions were Gaussian. ^ At the 
relatively fine resolution 4x4 array size, all seven pre- 
cipitation classes and the no-precipitation overcast data 
displayed bimodal distributions while at the relatively 
coarse 10 x 10 array size, all displayed an unimodal 
di stribution . 
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E. CONCLUSIONS 

The conclusions from this study of precipitation speci- 
fication are: 



• Varying the satellite data resolution from 434 nmis (4 
X 4 array size) to 2025 nmi2 (10 x 10 array size) 
results in a statistically significant difference in 
the representation of precipitation or no-prscipitat ion 
data. Variation in the distribution functions and the 
characteristic means and standard deviations with 
increasing cell size demonstrates systematic trends 
which may, with further study, provide improved basis 
for rain detection and/or quantification. 






For overcast cei 
tral threshold b 
tion is defined 
a cloud top tem 
and a cloud top 
ing reports. F 
a threshold base 
specifies a 50% 
the precicitatio 
lite resclution 



ling report 
as e d on a 5 
as extendin 
erature of 
emperature 
or overcast 
d on the al 
probability 
n probabili 
size . 



bi-spec- 

recipita- 



? 



s, a simple linear 
0% probability of prec; 
g from an albedo or 1.00 and 
290K to an albedo of 0.60 
of 210K for overcast ceil- 
and broken ceiling reports, 
bedo being greater than 0.80 
of precipitation. However, 
ties varied with the satel- 



• The precipitation and no -pr ecipit ation data sets can be 
differentiated by the mean of the means, the mean of 
the standard deviations, and the trend of the mean of 
the means. 



• Differentiation of convective precipitation from con- 
tinuous precipitation shows promise through the differ- 
ences in the infrared mean of the means and through the 
visual mean of the standard deviations. 



• Qualitative specification of light versus moderate/ 
heavy precipitation shows some promise through a 
threshold based on the relative shift toward colder 
cloud top temperatures and higher albedos of the moder- 
ate/heavy precipitation relative to the light precipi- 
tation. 



F. SUGGESTED FURTHER STUDY 



The recommendations for further study are: 



• The proposed data processing that would display the 
visual and infrared mean and standard deviation of the 
means (IV.C.1. ) should be investigated for both the 
overcast and the overcast and broken data to determine 
if the precipitation and no-precipitation data can be 
more distinctly separated. 
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• The appearance of bimodal distributions at the fine 
resolutions (6x6 and 4x4 array sizes) in rhe seven 
precipitation categories and in the no-precipitation 
overcast category should be invesrigar ed as a possible 
synoptic signature. 



• Formal discriminant analysis should be applied to the 
relevant statistics to yield confidence levels of indi- 
vidual visual and infrared satellite data pairs for 
precipitation specification. 



• The data should be further investigated at the 2 x 2 
array size to analyze this resolution in representing 
the data. 
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AEPENDIX A 



The Muench and Keegan (1979) normalization relates the nor- 
malized reflectivity, r^ , to the varying solar angle and 
maximum digital counts through the reflectance term, and 

the anisotropic scattering through the % term. Table XVI 
defines the symbols. Table XVII lists the geometric identity 
equations and Table XVIII list the normalization equations. 
Fig. 37 gives an example of the normalization applied to the 
stated location. 
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TABLE XVI 
List of Symbols 



Symbol 


Description 


Units 


c 


GOES video count number (0-63) 


dimensionless 


Co 


GOES video count number for perfect 
diffuse reflector and overhead sun 


dimensionless 


G 


Greenwich meridian time 


hours-minutes- 

seconds 


R 


Distance of earth to sun 


km 


R 

o 


Mean distance of earth to sun 


km 


d 


Julian date 


dimensionless 


h 


Hour angle 


radians 


r 


Cloud reflectivity 


dimensionless 


Y 


Arc-length observer to subsatellite 
point 


radians 


£ 


Declination of the sun 


radians 


C 


Zenith angle of the sun 


radians 


A 


Longitude 


radians 




Longitude of subsatellite point 


radians 


X 


Anisotropic scattering coefficient 


dimensionless 




Latitude 


radians 




Azimuth of the sun 


radians 


^2 : 


Azimuth of the satellite 


radians 
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TABLE XVII 

Basic Geometric Sat ellite- Earth Relationships * 

Declination : 

6 = 0.408 sin [(d~81) ♦ 2ir/365] 

Solar distance ratio: 

R/Rg s 1 _ 0.167 cos [(d-14) * 2ir/365] 

Hour angle: 

h = A + IT - G(hours) x ■ir/12 
Arc-length: 

cosy » cosCAg - A) cos<t 
Satellite azimuth: 

sin C'J'g - IT) = sin (Ag - A)/siny 
Solar azimuth angle: 

cose * sin4> sin6 cos<I> cos6 cosh 
Solar azimuth: 

sin4>^ * cos6 sinh/sinc 



Angles in radians 

TABLE XVIII 

Muench and Keegan Normalization Equations ♦ 

f =» • secc 

=« 1$2 - 

= COS^ cc? - 50) * 1.8) 

C 2 =■ 0.7 cos ((? - 22.5) * 4) * (1 - cos?) 

Cg = cos® ((A({> - 70) * 1.3) 

X * 1.0 + 0.05 * (1 + cos(2*0) + 0.20 * + C 2 ) * Cg 

» [1.09 - 2H1.09 - f * X *(R/Rq)^)/( 1 + cos^^^ 

* 

Angles in degrees 
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COUNTS 

10.0 15.0 20.0 25.0 30.0 35.0 40.0 45.0 50.0 55.0 




TIME (GMT) 



. Normalized Cloud Reflectivity as a Function 
Video Count and Time of Day for April 16 at 
and 740 W, with Calibration Count. (C ) equal 
60. (Example taken from Musnch and°Keagan, 
1979.) 
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Figure 37 
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